Crossed Andreev reflection in superconducting graphene spin- valves: Spin-switch effect 
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We consider the non-local quantum transport properties of a graphene superconducting spin-valve. It is shown 
that one may create a spin-switch effect between perfect elastic co-tunneling (CT) and perfect crossed Andreev- 
reflection (CAR) for all bias voltages in the low-energy regime by reversing the magnetization direction in 
one of the ferromagnetic layers. This opportunity arises due the possibility of tuning the local Fermi-level 
in graphene to values equivalent to a weak, magnetic exchange splitting, thus reducing the Fermi surface for 
minority spins to a single point and rendering graphene to be half-metallic. Such an effect is not attainable in 
a conventional metallic spin-valve setup, where the contributions from CT and CAR tend to cancel each other 
and noise-measurements are necessary to distinguish these processes. 

PACS numbers: 74.25.Fy,74.45.+c,74.50.+r,74.62.-c 
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I. INTRODUCTION 

Quantum entanglement' describes a scenario where the 
quantum states of two objects separated in space are strongly 
correlated. These correlations can be exploited in emerging 
technologies such as quantum computing, should one be able 
to spatially separate the entangled objects without destroy- 
ing the correlations. In a broader context, quantum entangle- 
ment could prove to be of practical importance in the fields 
of spintronics^ and information cryptography^ It also holds a 
considerable interest from a purely fundamental physics point 
of view, prompting some of the more philosophically inclined 
discussions related to quantum theory and causality. 

Superconductors have been proposed as natural sources 
for entangled electrons'*'^, as Cooper pairs consist of two 
electrons that are both spin and momentum-entangled. The 
Cooper pair can be spatially deformed by means of the crossed 
Andreev reflection (CAR) process in superconducting het- 
erostructures. In this scenario, an electron and hole excitation 
are two separate metallic leads are coupled by means of An- 
dreev scattering processes at two spatially distinct interfaces. 
Unfortunately, the signatures of CAR are often completely 
masked by a competing process known as elastic co-tunneling 
(CT) which occur in the same type of heterostructures. In fact, 
the conductances stemming from CT and CAR may cancel 
each other completely^, thus necessitating the usage of noise- 
measurements to find fingerprints of the CAR process in such 
superconducting heterostructures. 

Recently, graphene^ has been studied as a possible arena for 
CAR-processes. In Ref.**, it was shown how a three-terminal 
graphene sheet containing n-doped, p-doped, and supercon- 
ducting regions could be constructed to produce perfect CAR 
for one particular resonant bias voltage. Also, the signatures 
of the CAR process in the noise-correlations of a similar de- 
vice were studied in Ref.''. However, the role played by the 
spin degree of freedom in graphene devices probing non-Iocal 
transport has not been addressed so far. This is a crucial point 
since it might be possible to manipulate the spin-properties of 
the system to interact with the spin-singlet symmetry of the 
Cooper pair in a fashion favoring CAR. 
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FIG. 1: (color online) Proposed experimental setup for the spin 
switch effect between crossed Andreev reflection and elastic co- 
tunneling. Ferromagnetism and superconductivity are induced by the 
proximity effect to a host material. The induced exchange fields in 
the non-superconducting graphene regions are oriented either paral- 
lel or antiparallel with respect to each other. In the parallell align- 
ment, the density of states vanishes for both normal Andreev reflec- 
tion and crossed Andreev reflection processes, such that only elastic- 
cotunneling contributes to non-local transport. In the anti-parallel 
alignment, the density of states vanishes for both normal Andreev 
reflection and elastic co-tunneling, leaving only crossed Andreev re- 
flection as the non-local transport channel. 



In this paper, we show that precisely such an opportunity 
exists - it is possible to obtain a spin-switch effect between 
virtually perfect CAR and perfect CT in a superconducting 
graphene spin valve. In contrast to Ref.^, this effect is seen 
for all bias voltages in the low-energy regime rather than just 
at one particular applied voltage difference. The key obser- 
vation is that the possibility of tuning the local Fermi-level 



2 



to values equivalent to a weak, magnetic exchange splitting 
in graphene renders both the usual Andreev reflection process 
and CT impossible. In contrast, this opportunity does not ex- 
ist in conventional conductors where the Fermi energy is large 
and of order ©(eV). We show that graphene spin valves pro- 
vide a possibility for a unique combination of non-local An- 
dreev reflection and spin-dependent Klein tunneling Our 
model is shown in Fig. [T] where ferromagnetism and super- 
conductivity are assumed to be induced by means of the prox- 
imity effect" '^ to leads with the desired properties. A similar 
setup was considered in Ref.'^*, where the magnetoresistance 
of the system was studied. 

We organize this work as follows. In Sec. [Ill we establish 
the theoretical framework which will be used to obtain the 
results. In Sec. [Till we present our main findings for the non- 
local conductance in the graphene superconducting spin- valve 
with a belonging discussion of them. Finally, we summarize 
in SecHVl 



II. THEORY 

We consider a ballistic, two-dimensional graphene structure 
as shown in Fig. [T] In the left ferromagnetic region a; < 0, 
the exchange field is h = /iqz, while it is h = ±hoz in the 
right ferromagnetic region x > L. In the superconducting 
region Q < x < L, the order parameter is taken to be con- 
stant with a real gauge A = Aq. To proceed analytically, we 
make the usual approximation of a step-function behavior at 
the interfaces for all energy scales, i.e. the chemical poten- 
tials {iJ,F,fJ's}, the exchange field ho, and superconducting 
gap Aq. This assumption is expected to be good when there 
is a substantial Fermi-vector mismatch between the F and S 
regions, as in the present case. To make contact with the ex- 
perimentally relevant situation, we assume a heavily doped S 
region satisfying ^5 ;» 

We use the Dirac-Bogoliubov de Gennes equations first em- 
ployed in Ref.'"^. For quasiparticles with spin a, one obtains 
in an F|S graphene junction: ' ^I'^i'^i^^''^ 



and in this case reads 
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and . " . denotes a 2 x 2 matrix. Here, we have made use 
of the valley degeneracy and p is the momentum vector in 
the graphene plane while a_ is the vector of Pauli matrices 
in the pseudospin space representing the two A, B sublat- 
tices of graphene hexagonal structure. The superconducting 
order parameter A(x) couples electron- and hole-excitations 
in the two valleys (=f) located at the two inequivalent corners 
of the hexagonal Brillouin zone. The spinor describes the 
electron-like part of the total wavefunction 



(4) 



while u^'^ = Tu"^ . Here, ^ denotes the transpose while T is 
the time-reversal operator. 

From Eq. ([T]), one may now construct the quasiparticle 
wavefunctions that participate in the scattering processes^. 
We consider positive excitation energies £ > with in- 
coming electrons of n-type, i.e. from the conduction band 
e = wfIpI — (we set — \ from now on). The incom- 
ing electron from the left ferromagnet may either be reflected 
normally or Andreev-reflection (AR). In the latter process, it 
tunnels into the superconductor with another electron situated 
at (— e), leaving behind a hole excitation with energy e. The 
scattering coefficients for these two processes are and Vh, 
respectively, and the total wavefunction may thus be written 
as: 
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where we have defined the wavevectors 



(5) 



(6) 



We have omitted a common factor e'^"" for all wavefunctions. 
Similarly, assuming that the charge carriers in the right ferro- 
magnetic region are also of the n-type, we obtain: 




ip^^ cos X 



Q-ipf cose^'x^ 



(7) 



It should be noted that the AR hole is generated in the con- 
duction band \f e — jip ~ crh^ > (retro- AR), whereas it is 
generated in the valence band otherwise (specular-AR). The 
± sign above refers to parallell/antiparallell (P/AP) magneti- 
zation configuration. 

We assume that the superconducting region is heavily 
doped, /i5 ^ /iF + ^0. which causes the propagating quasi- 
particles to travel along the x-axis since the scattering angle in 
the superconductor satisfies ^5 0. We obtain the following 
wavefunction (A = ±1): 
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where k = -^/Aq — while 

P - acos(£/Ao) (9) 

for subgap energies |e| < Aq and 

f3 = -iacosh(e/Ao) (10) 

for supergap energies \e\ > Aq . 

It is important to consider carefully the scattering angles in 
the problem. Since we assume translational invariance in the 
y-direction, the y-component of the momentum is conserved. 
This gives us 

sin 61 = sin 61^ = pf" sin 9%". (11) 

It is clear that the angle of transmission for the electrons in 
the right ferromagnet is equal to the angle of incidence when 
the magnetizations are P, i.e. 6'^ — 9. Also, one infers that 
there exists a critical angle above which the scattered waves 
become evanescent, i.e. decaying exponentially. This may 
be seen by observing that the scattering angles exceed 7r/2 
(thus becoming imaginary) above a certain angle of incidence 
9. For instance, the AR wave in the left ferromagnetic region 
becomes evanescent for angles of incidence 9 > 9'^^, where 
the critical angle 9 = 9'^g^ is obtained by setting 0^ = tt/2 in 
the equation 

p^sin6' =p^sin6'^, (12) 

expressing conservation of momentum perpendicular to the 
interface. One finds that: 

= |asin[(£ - fip + (yha)/{e + ^ip + (jho)]\. (13) 

Thus, AR waves in the regime 9 > \9^\ do not contribute 
to any transport of charge. A similar argument can be made 
for the transmitted electron wave-function in the right ferro- 
magnetic region, corresponding to the CT process, where the 
critical angle for this process becomes 

91j = |asin[(e + /iF ±c^/io)/(e + MF + ct/io)]|- (14) 

In the P configuration, the CT process thus always contributes 
to the transport of charge. Finally, the contribution to transport 
of charge from CAR comes from the hole-wave function in 
the right ferromagnetic region, which becomes evanescent for 
angles of incidence above the critical angle 

^'cAR > |asin[(e - a^f ± cr/io)/(e + + ct^o)]|- (15) 

In the P configuration, this criteria is the same as the vanishing 
of local AR expressed by Eq. (fT3l) . 

III. RESULTS AND DISCUSSION 

Intuitively, one might expect that the most interesting phe- 
nomena occur when the exchange field Kq is comparable in 
magnitude to the chemical potential ^ip. If ^ip » h^, the ef- 
fect of the exchange field should be minor and the AR is never 
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FIG. 2: (color online) Plot of the conductance for CT processes 
Gct/Gf versus bias voltage in the upper panel and versus length 
of the S region in the lower panel. Here, we consider the P alignment 
and jj,F ~ ho such that Gcar — *■ 0. 



specular. In contrast, the situation becomes quite fascinating 
when we consider the case fip = ho under the assumption 
of a doped situation ^p ^ {e, Aq). First of all, the incom- 
ing quasiparticles from the left ferromagnetic region are com- 
pletely dominated by the majority spin carriers a ^1, since 
the density of states (DOS) for a =J, electrons vanishes at the 
Fermi level. Since ^J,p — Iiq, the AR process is suppressed for 
all incoming waves as 0^,^ 0. We now show how the fate 
of the cross-conductance in the right ferromagnetic region de- 
pends crucially on whether the magnetization configuration is 
P or AP. In the P configuration, we see that 0, which 

means that the transport is purely governed by the CT process. 
In the AP configuration, we see that 9^j 0, which means 
that the transport is mediated purely by the CAR process. This 
suggests a remarkable spin-switch effect - by reversing the 
direction of the field in the right ferromagnet, one obtains an 
abrupt change from pure CT to pure CAR processes mediat- 
ing the transport of charge. In each case, there is no local AR 
in the left ferromagnetic region. In the standard metallic case, 
the distinct signatures for the CT and CAR contributions are 
masked by each other, and it becomes necessary to resort to 
noise-measurements in order to say something about the con- 
tribution from each process. In the present scenario, we have 
showed how it is possible to separate the two contributions 
directly by a simple spin-switch effect which is commonly 
employed in experimental work on F|S heterostructures. 

Let us now evaluate the conductance in the P and AP con- 
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figuration quantitatively by using 



Gcar/Gf = V(g"/Gf) / decose\th\^, m 

J-7r/2 



where we have introduced 



(17) 



as the spin-(7 normal-state conductance that takes into account 
the valley degeneracy, in addition to 



Gp — G^ -\- G . 



The density of states is determined by 



N''{e) = \e + ^iF + ah\W/{-KVY), 



(18) 



(19) 



where W is the width of the junction. The expression for 
GcT is obtained by replacing th with in Eq. (fT6b . Since 
we here consider the case fip = fiQ and ho 3> (e, Aq), 
the formulas for the Gcar and Gct may be simplified since 
G_ ^ G+. Also, since the DOS vanishes for minority spins 
for the injected electrons, only cr =1 contributes for incoming 
electrons. The crucial point here is that in the P alignment, 
Gcar and Gct such that 



= 1, 



(20) 



while in the AP alignment Gcar and Gct such that 



\thr = 1- 



(21) 



In the actual numerical calculations, we use h^/Ao = 50 and 
lig/Af) — 500. Assuming a value of Aq = 0.1 meV for 
the proximity-induced gap, this corresponds to an exchange 
splitting of ft-o = 5 meV in the F regions and a doping level 
/is = 50 meV in the S region, which should be experimen- 
tally feasible^ and well within the range of the validity for 
the linear dispersion relation in graphene. In Fig. |2] we plot 
the cross-conductance Gct / Gp in the P alignment both as a 
function of bias voltage and width of the S region. The same 
thing is done for Gcar/Gf in the AP alignment in Fig. [3] In 
both cases, the magnitude of the conductance varies strongly 
when considering different widths L due to the fast oscilla- 
tions which pertain to the formation of resonant transmission 
levels inside the superconductor. Also, it is seen that while the 
CT process is favored for short junctions L/£^ ^ 1, the CAR 
process is suppressed in this regime in favor of normal reflec- 
tion. Upon increasing the junction width, the CT conductance 
drops while the CAR conductance peaks at widths L ^ ^. 
The remarkable aspect is that it is possible to switch between 
these two scenarios of exclusive CT and exclusive CAR sim- 
ply by reversing the direction of magnetization in one of the 
ferromagnetic layers. 

In order to obtain analytical results, we have assumed that 
the Coulomb interaction and charge inhomogeneities may be 
neglected. It would be challenging to obtain a truly homoge- 
neous chemical potential in a graphene sheet, and electron- 
hole puddles appear to be an intrinsic feature of graphene 
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FIG. 3: (color online) Plot of the conductance for CAR processes 
Gcar/Gf versus bias voltage in the upper panel and versus length of 
the S region in the lower panel. Here, we consider the AP alignment 
and jj,F = ho such that Gct 0. 



sheets^^. Moreover, it has been speculated that such charge 
inhomogeneities may play an important role with regard to 
limiting the transport characteristics of graphene^^ near the 
Dirac points. However, for our purposes this is actually bene- 
ficial - it is precisely the suppression of charge and spin trans- 
port at Fermi level for the Andreev reflection and co-tunneling 
process which renders possible the spin-switch effect. There- 
fore, we do not expect that the inclusion of charge inhomo- 
geneities should alter our results qualitatively. Finally, we 
note that since the spin of the charge-carriers in each of the 
non-superconducting graphene sheets are practically speaking 
fixed due to the vanishing DOS for minority spins, the spin- 
switch effect for CAR and EC predicted in this paper can not 
be directly related to entanglement. Nevertheless, it consti- 
tutes a clear non-local signal for quantum transport which can 
be probed experimentally, and should be helpful in identifying 
clear-signatures of the mesoscopic CAR phenomenon. 



IV. SUMMARY 

To summarize, we have considered non-local quantum 
transport in a graphene superconducting spin-valve. We have 
shown how one may create a spin-switch effect between 
perfect elastic co-tunneling and perfect crossed Andreev- 
reflection for all applied bias voltages by reversing the magne- 
tization direction in one of the ferromagnetic layers. The basic 
mechanism behind this effect is that the local Fermi-level in 
graphene may be tuned so that the Fermi surface for minority 
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spins reduces to a single point in the presence of a weak, mag- and 167498A'^30 (STORFORSK). 
netic exchange spHtting. This is very distinct from the equiv- 
alent spin valve structures in conventional metallic systems, 
where noise-measurements are required to clearly distinguish 
between these processes. 
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